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Abstract 



In the many-body theory of particles with long range interaction potential, 
such as the potential of Newtonian gravity and its relativistic generalization, 
Newton-Birkhoff theorem plays a very important role. This theorem is vi- 
olated in the case of short range interaction potential, such as the Yukawa 
potential in nuclear physics. We discuss the relevance of this in the many- 
body treatment of heavy nuclei. In particular, using techniques similar to 
those in Newtonian gravity, we calculate analytically the mean-field potential 
in heavy nuclei. The mean-field potential thus calculated depends on well 
known quantities such as inverse pion mass, pion-nuleon coupling constant, 
nuclear radius and density. The qualitative features of the mean-field poten- 
tial thus obtained is similiar to the well known phenomenological Wood-Saxon 
potential. 
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In this paper, using simple analytical calculations, we intend to bring out some fea- 
tures of motion of nucleons in a heavy nucleus [1,2]. Our calculations are based on two 
phenomenologically well known facts: The density of nucleons in heavy nuclei is almost con- 
stant and the inter-particle interactions in nuclei is given by the Yukawa potential. These 
assumptions seem to be sufficient for the qualitative understanding of the motion of nucle- 
ons in heavy nuclei. The central ingredient in our calculations is that the Yukawa potential 
violates Ncwton-Birkhoof theorem. In order to explain the role played by the violation of 
this theorem by the Yukawa potential in nuclear physics, we first formulate and proof the 
theorem in Newtonian gravity. The proof reveals why such a theorem can not be valid for 
Yukawa potential. In the many-body treatment of heavy nuclei with inter-particle Yukawa 
interactions, we consider the nucleon density to be constant. In order to justify this as- 
sumption, in the following section, we argue why a heavy nucleus can be treated as system 
of constant density. We show that the constant density has more to do with the Pauli's 
exclusion principle and short range repulsive core of the nucleons rather than any ergodicity 
or thermodynamic equillibrium [3-7]. In the next two sections, we derive the mean field 
potential and indicate how the surface tension of heavy nuclei emerges from the violation 
of Newton-Birkhoff theorem. The mean field potential thus obtained qualitatively agrees 
with the phenomenological Wood-Saxon potential [8]. At the end, we make some concluding 
remarks. 



I. THE NEWTON-BIRKHOFF THEOREM 

Newton-Birkhoff theorem in the theory of gravitation is a very important theorem. The 
violation of this theorem by Yukawa potential is the central ingredient in our calculation of 
the surface tension and the mean field of heavy nuclei. Therefore, in this paper we provide 
a discussion of this theorem in Newtonian gravity in detail. Let us consider a spherically 
symmetric ball of mass M, mass density p and radius R. 

Theorem (Part-I): The potential due to the spherical ball at a distance r, where r > R is 
the same as the potential due to point particle of mass M placed at the centre of spherical 
ball. 

Let us consider a spherical shell of outer radius R and inner radius b, having constant mass 
density p. 

Theorem (Part-II): The gravitational potential due the spherical shell at any point r, where 

r < 6 is constant. Therefore, the force at the point r is zero. 

This theorem also holds in general theory of relativity and is called the Newton-Birkhoff 
theorem. Bellow we provide proof of the theorem in the contex of Newtonian Gravity. The 
proof reveals why the theorem is so special to Newtonian gravitatinal potential and does 
not hold for other potentials. 

We want to calculate the gravitational potential and potential energy of a spherically 
symmetric ball of radius R and constant density p.The two-body Newtonian gravitational 
potential energy is given by, 

U,^-^^ (1) 
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where r and r' are the locations of the point masses mi and m2 respectively. The Newtonian 
gravitational potential at a point r due to a point mass m at point r' is given by 

^ , . Gui , „ , 

*iv(r) = (2) 
|r — r I 

For a mass distribution with density p(r'), the formulas above take the following forms, 



*„(f, = _G/^(^ (3) 



In spherical co-ordinate system the formula for $Ar(r) gives, 

/ / /• r''^dr'sineded(f) 
*^^^) = -^^yyy (r^ + r--2rrW)V2 

In the expression above, the integral over ranges from to 27r and for 9, it ranges from 
to TT. The integral over gives a factor of 27r. The integral over 9 can be carried out by 
making the following change of variable 

^2 ^ + r'^ - 2rr'cos9- ^d^ = rr'sin9d9 ; 

The integration over ^ ranges from C — (k ~ ^1) to ^ — {r + r') After carrying out 
integrations over and 9, we obtain 

$iv(r) = - 27rGp y" [(r + r') - \r - r'\]'^dr' (6) 
For r > a (> r'), and < r' < a 

$jv(r) = -2'kGp [(r + r')-(r-r')]-dr' ^-477— r'^dr' ^ — ^ = (7) 

JQ r r Jo or r 

For r < 6 (< r'), and b < r' < R 

$iv(r) = - 27rGp /"[(r + r') - (r' - r)]-dr' = -AnGp r'dr' = -27rGp(i?=^ - 6^) (8) 
Jo r Jb 

which is constant independent of r . Note that, in the equation above, the cancellation of 
r from the expression is due to the special form of the Newtonian potential. Since force 
is given by the gradient of the potential with respect to r, F = — Vf?$Af(r) = 0, therefore 
there is no gravitational force within the spherical shell. This demostrates that in a classical 
theory this constant potential is dynamically irrelevant. However, it should be noted that 
in a quantum theory this "constant potential well" is relevant and its' physical effects can 
be observed. 
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II. CONSTANT DENSITY OF NUCLEONS IN HEAVY NUCLEI 



In this section we discuss why the density of nucleons in heavy nuclei can be taken to 
be constant. We do not provide here a formal proof of our claim. We rather argue that it 
follows very natuarally from the short range character of the attractive potential, Pauli's 
exclusion principle and the repulsive hard core potential of nucleons. Our claim is essentially 
based on self consistency arguments. However, it is important to note that we do not require 
any ergodicity or thermodynamic equillibrium for the constancy of the density of nucleons 
in the nuclei. 

In order to maintain simplicity in our exposition and to avoid complications, we consider 
only one type of micleon and assume that it does not carry any electric charge. The attactive 
interaction potential between a pair of nucleons is given by the short range Yukawa Potential, 

V = -g (9) 

r 

2 

where g is the nuclear charge ( ^ is the dimensionless nuclear fine structure constant ) and 
jjL characterizes the range of the attractive potential. We assume that d is the radius of the 
core of the hard core potential. Wc assume that, i?, is the radius of the nucleus and is 
the number nucleons in it. At first, we consider the nuclear matter limit of the nucleus. 
This means that both R and are large. Yukawa potential is a short range potential, and 
therefore, each nucleon interacts with nucleons in a small neighbourhood around itself. The 
number of nucleons in its neighbourhood is proportional to the density of nucleons Up . Let 
the constant of proportionality be / . Note that / is the volume charaterising the range of 
the potential and can be taken to be / ^ Let us represent by —a the potential energy 
of attractive interaction between two nucleons. Therefore, the potential energy of attractive 
interaction of a single nucleon with its neighbours is —f ripa , and of A^ nucleons with their 
neighbours is —fNupU. In the next section, we will show that fa — Ang'^/iJ,'^ (see the end 
of section-Ill). 

Wc assume that the kinetic energy of the nucleons in the nucleus is solely due to the zero 
point energy. Therefore, using Pauli's exclusion principle we obtain, K.E. ~ ^ x "^^r-- 
Taking Up ~ 47r/3 r3 ) we write the total energy E[R) as function of R as, 

E(R) pa — X — — - fNrioa = — x — V (10) 

We look for the variational (QM) ground state of the system. This can be done by finding 
the minimum of the function E(R). By imposing the condition, ^ = 0, wc obtain all the 
cxtrcma. It turns that the function has a minimum at i? = (£"(-/? = 0) = — cxo) and a 
maximum when the interparticle distance (= is a few fermi. The function goes to zero 
at infinity. Therefore, there is no minima for finite value of R. Such a system would either 
fiy apart or collapse to a point. The ground state of the system is unstable. It is here that 
the hard core potential comes to play an important role. Because of the hard core potential 
the collapse would stop when = d . Therefore, 

" 2m R^ AnR^ ~ 2m d^ 47r# ^^^^ 
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and the per particle binding energy 



3fa 



(12) 



47r# 



This leads to the emergence of constant density. This picture resembles very much like those 
in packed granular material. However, our arguments are exact only in the limit of nuclear 
matter. For finite nuclei, however heavy it may be, volume occupied by the nuclcons will be 
a bit larger than the hard core radii. Pauli's exclusion principle adding a soft core in addition 
to the already existing hard core. This soft core is not a scattering potential. Because of 
Pauli's exclusion principle , the nucleons rather avoid each other. It is a granular material 
in which the particle fail to touch other by the virtue of Pauli's exclusion principle. But the 
gap is limited and this does not allow the particles to pile up at any one single place thereby 
mainintaining a constant density profile. 

III. DERIVATION OF THE MEAN FIELD POTENTIAL IN HEAVY NUCLEI 

We consider heavy spherical nucleus consisting of N nucleons and of radius R . The 
density of nucleons in the nucleus, Up, is considered to be constant. Since the system is 
spherically symmetric, the potential at any point on a layer at r will be the same. In general 
the potential on any layer at r will consist of two parts- the potential due to the spherical 
ball of radius r — e and the potential due to the spherical shell of inner radius r + e and 
outer radius R. We will calculate these two contributions separately and then add them up. 
This amounts to calculating the potential, out side a spherical ball and the potential , 
$2, inside a spherical shell exactly in the way we carried out the calculations in Newtonian 
gravity. Let us write the general expression for the potential, 



In the expression above, the integral over ranges from to 27r and for 6, it ranges from 
to IT. The integral over gives a factor of 27r. The integral over 6 can be carried out by 
making the following change of variable 



The integration over ^ ranges from C = ~ '"'I) to ^ — {r + r') After carrying out 
integrations over and 9, we obtain 




(13) 



= + r'^ — 2rr'cos9; ^d^ = rr'sinOdO ; 




(14) 



For r > a {> r'), and < r' < a 




(15) 



After carrying out the integration, we obtain 
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4:7rgnp e , / \ 1 , / m 

<I>i(r, a; r > a) = [a cosh[(ia) sinh{na)\ (16) 

/J, r /J, 

This is the potential out side the spherical ball of radius a at a distance r from its' centre. 
Let 

AiTTQTl 1 

A(a) = coshiixa) sinhiiia)] 

1^ // 

This allows us to write the potential outside the spherical ball as 



e 



$i(r, a; r > a) = -A{a) 

which has the same form as the Yukawa potential except for the fact that, now, the nuclear 
charge is A{a) and not g.To obtain the same potential as that of the spherical ball at distance 
r, r > a we need to place at centre of the ball a point nucleus of nuclear charge equal 
to A{a). We can estimate the value of A{a) by writing the series expansion of hyperbolic 
cosine and sine functions in the expression of A{a). This gives, 

r 1 1 °° 277 1 

A{a) = 7V^[l + -(H' + ^(/^«)' + E 3 X (17) 

where N — is the number of nucleons inside the spherical ball of radius a. When 

radius a is large, the quantity A is much much larger than the total nucleon charge Ng 
inside the spherical ball. Therefore, to obtain the same potential as that of the spherical 
ball containing the amount of nuclear charge Ng ( at distance r, r > a ), we need to place 
at centre of the ball a point nucleus of nuclear charge equal to A which is larger than Ng. 
This clearly demostrates that the part-I of the Newton-Birkhoff theorem is violated by the 
Yukawa potential. However, in the limit of — > 0, we recover result similar to the case of 
gravity. 

For r < 6 (< r'), and h <r' < R 

27rgn '"^ 
/ir 

After carrying out the integration, we obtain 



$2 (r, 6 ; r < 6) = -^I^ [ dr' r'fe-'^^^'-^) - e~^^'+'-'^] (18) 

Ltr Jb 



jjL'^ r jji jjL 

$2 {r,b -jT < b) is the potential inside a hollow spherical shell, F = — V,^$2 (r, 6 ; r < 6) 7^ 0, 
and therefore, there is force inside a hollow spherical shell. This violates part-II of the 
Newton-Birkhoff theorem. However, in the limit of ^ 0, we recover result similar to the 
case of gravity. 

It is now easy to calculate the potential at any r within the nucleus. Assuming, 

B(x)^^^(x+-)e-^'' (20) 
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We can write the total potential at r, a < r < 6 as 



-A{a) B{b) 



sinh{fir) 



+ B{R) 



sinh{fir) 



(21) 



When the radius of the inner sphere, a , approaches r from bellow and the inner radius 
of the shell, b , approaches r from above , the first two terms of the equation above add 
up to a constant independent of r , the third term remains unchanged and we obtain 



1^' 



1-(R+-) e 



sinh{iir) 



This is the mean field potential of the nuclcons at any point r inside the nucleus. The 
mean potential outside the nucleus is given by 



Angnp e 



fir 



R cosh(iiR) sinhi^R) 



(22) 



The potential energy of a single nucleon in the mean mean field inside and outside the 
nucleus is given by the following equations. 



U{r < R) 



i-(r + 



1 X sinhjur) 
ji' r 



(23) 



U{r > R) 



R cosh{iiR) sinh{iiR) 

A* 



(24) 



In the formula above the values of the the physical constants arc the following: pion- 
nucleon coupling constant, g^/ch = 0.081 , conversion constant, ch = 197.328, range of 
two-body Yukawa potential, l/jj, = 1.4094 fm, radius of nucleus of atomic number A, 
R — 1.07 A^/^ fm, density of nucleons in nucleus, Up = 0.0165//m^. We provide plots of 
the potential for three different nuclei in Fig.(l). The quahtative feature of this potential 
is similiar to the phenomenological Wood-Saxon potential [8]. The qualitative features is 
similiar even near the boundary of the nucleus. This is a bit of a surprising result because 
the nucleon density near the boudary of the nucleus decreases drastically where as in our 
derivation we have assumed constant density through out the nucleus. This suggests that 
the mean field potential is not very sensitive to the distribution of nucleons near the bound- 
ary. 

From eq.(23), it is now easy to infer that in the nuclear matter limit, i.e., when R is very 



large, U (r) 
47r^7/x2 . 



4:Trg Up 



Therefore, the quantity, fa introduced in section-II is equal to 



IV. THE ORIGIN OF SURFACE TENSION IN HEAVY NUCLEI 

In this section we sketch some arguments to show that the surface tension in heavy nuclei 
has its' orgin in the violation of Newton-Birkhoff theorem. Actual calculation of the surface 
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tension will require the true density profile of the nucleons in the nucleus. However, in the 
frame work of the present paper the density is assumed to be constant, and therefore, we 
will not attempt to calculate the surface tension. We rather point out that , even under 
the assumption of constant density, one can observe the origin of surface tension. The true 
density profile infiuences only the numerical value of the surface tension. 

To make the role played by the violation Newton-Birkhoff theorem more transparent, 
we briefly recall some of the derivations. The total potential at any point, r, in the region 
between a sphere of radius a and a spherical shell of inner radius b and outer radius R 
( a < r < 6 ) is given by 

^r) = -^(a)^ - B{bf-^^^ + B{R) 

r r r 

When the radius of the inner sphere, a , approaches r from bellow and the inner radius 
of the shell, h , approaches r from above , the first two terms of the equation above add 
up to a constant independent of r , the third term remains unchanged and we obtain 



Mr) = -^[l-(i?+-)e 



Force F , on a single nucleon at any point r , r < R is given by. 



1^ ' 1^ 



and for r , r > R 



F = -gV^outir) = + r ^ ^^^^^^^m _ Isinh(iiR) In, (26) 

where Hr is a unit vector directed along the radius towards the centre of the nucleus. 
The depth of the mean field potential varies very slowly in the central region (r ?s 0) of the 
nucleus ( see Fig.(l) ). However, near the outer surface, there is a sharp decline in the depth 
of the total potential. Therefore, the force will have very small in magnitude (practically 
zero) in the central region but a broad peak near the bounday of the nucleus. This is what 
we observe in Fig. (2), where we have plotted the total force, F as a function of r. The 
direction of this force at any point r near the boundary is towards the centre along the 
radius connecting the point. Every nucleon in the outer shell is, therefore, radially pulled 
towards the centre of the nucleus. Since the force is practically zero in the central part of 
the of the nucleus, there is a net pressure difference between the outer shell and the inner 
core. This is the origin of surface tension in nuclei. We beheve that the surface tension in 
most normal liquids have similar origin . Our believe is based on the fact that potentials 
of the type — 1/r" which describe the properties liquid state very well, can be obtained by 
the exchange of continuum of Yukawa modes [9,10] with some spectral weight p(/i). In fact, 
the most commonly used Leonard- Jones potential in liquids can be obtained with spectral 
weight p{ii) = i/". 

1 1 /-oo p-l^r 

= / ^/^/^ (27) 

r° 4! 7o r 
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V. CONCLUDING REMARKS 



In this paper, we have tried to emphasise the importance of the violation of Newton- 
Birkhoff theorem by Yukawa interaction for the emergence of mean field shell model potential 
as well as the surface tension in heavy nuclei. The only input from the nuclear physics that 
we have used is the constant density profile of nucleons in heavy nuclei. Some very basic 
and general assumptions, such as Pauli's exclusion principle, short range character of the 
attractive Yukawa potential and the existence of a repulsive hard core in nucleons, ensures the 
emergence of constant density profile in heavy nuclei. We have derived the mean field shell 
model potential analytically. This potential qualitatively agrees wtih the phenomenological 
Wood-Saxon potential. Our formula explicitly shows how the mean field potential depends 
on pion-nucleon coupling constant, range of the two-body Yukawa potential, nuclear size and 
the density. Solution of the Schroedinger equation with the mean field potential so derived 
leads to the emergence of shell model structure. We have also pointed out how the violation 
of the Newton-Birkhoff theorem leads to the emergence of surface tension in heavy nuclei. 
Our derivations, although semi-quantitative in nature, strongly suggest that the qualitative 
features which are attributed either to shell model or liquid drop model can be derived from 
the Yukawa interaction of nucleons having repulsive hard core, and some basic principles of 
quantum mechanics. 
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Figure 1 : Plot of potential energy of a single nucleon in the mean-field, 
inside and outside the nucleus, as a function of radial distance from the 
center of the nucleus. 
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Figure 2: Plot of mean-field force F on a single nucleon inside and outside 
the nucleus as a function of radial distance /-from the centre of the nucleus. 

The minima in the plot occur at the radius R=1 .07 A^'^ for the respective nuclei. 



